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The propagation of an initially localized perturbation via an interacting many-particle Hamilto- 
nian dynamics is investigated. We argue that the propagation of the perturbation can be captured 
by the use of a continuous-time random walk where a single particle is traveling through an active, 
fluctuating medium. Employing two archetype ergodic many-particle systems, namely (i) a hard- 
point gas composed of two unequal masses and (ii) a Fermi-Pasta-Ulam chain we demonstrate that 
the corresponding perturbation profiles coincide with the diffusion profiles of the single-particle Levy 
walk approach. The parameters of the random walk can be related through elementary algebraic 
expressions to the physical parameters of the corresponding test many-body systems. 
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The transport properties of many-particle systems are 
of salient interest in diverse contexts, ranging from foun- 
dations of thermodynamics to the transduction of infor- 
mation on the nanoscale. The collective evolution of 
^ 1 interacting particles creates a dynamical "tis- 
sue," whose properties depend not only on the Hamilto- 
nian of the system, but also on the state of the system 
itself. An objective of primary interest is how the sys- 
tem responds to the perturbation that locally affects its 
dynamics. The answer then provides direct insight into 
collective energy, correlation, and information transport 
in extended nondissipative media [ll-Q]- 

Consider the situation of a many-particle system at mi- 
crocanonical equilibrium, when at the initial time t ~ 
one of the particles receives some external local pertur- 
bation. The system gains a small amount of perturba- 
tion energy, which is conserved due to Hamiltonian evo- 
lution of the system. However, the perturbation does 
spread as the perturbation energy is shared by a con- 
stantly growing number of particles. One of the main 
features of the spreading is the finite velocity of per- 
turbation propagation, < oo 0, i, 0. Therefore, 
similar to relativistic diffusion theories [8[, an effective 
"light cone" occurs Q, such that at the given time t the 
perturbation is almost completely confined to the inter- 
val [—vot,vot\ The fundamental fact of the cone's 
existence assumes the characteristics of a mathematical 
existence theorem [ll|; its strength therefore is the 
generality when dealing with many-particle systems in 
rather arbitrary situations. A pronounced weakness of 
this mathematical approach, however, becomes evident 
whenever one attempts to implement the theory for a 
particular system to obtain qualitative results on an an- 
alytical level of description. Therefore, more applied ap- 
proaches, like the diffusion formalism, might become of 
beneficial use. Indeed it is appealing to consider the said 
perturbation spreading as a certain (not yet known) one- 
dimensional diffusion process, and to quantify it with the 
mean square displacement and other related attributes 




FIG. 1: (color online) The standard continuous-time random 
walk (above the slab) and random walk through an active 
medium (in the slab). The first walker performs flights with 
constant velocity ±vo, while the second one is subjected to 
constant scattering, so that its velocity fluctuates. The dura- 
tion of single flight, r„, is a random variable drawn from the 
probability density function ipiT). 



[H [B-Ql- As it is known, however, it is impossible to 
describe this process by any known macroscopic, norm- 
preserving normal diffusion equation. Conventional dif- 
fusion equations knowingly lead to infinite propagation 
speeds [8| and therefore are incompatible with the exis- 
tence of a causal cone. Thus the perturbation kinetics 
should be ultimately considered on the microscopic level 
corresponding to the random walk approach. 

With this Letter we employ the microscopic single- 
particle random walk process in order to evaluate the 
evolution of perturbations in one-dimensional, ergodic 
many-particle systems. In doing so we address the two 
challenges: (i) How is the perturbation distributed within 
the cone, and (ii) what are the shapes of the cone fronts? 
By using two renowned many-particle chains we demon- 
strate that the random walk model accurately describes 
the perturbation spreading in these interacting many- 
body systems. We show that the walker-media interac- 
tions are responsible for the observed shape of the causal 
cone and predict the universal scalings for the perturba- 
tion profile and its corresponding fronts. 

Model setup. The continuous-time random walk 
(CTRW) formalism [l^l has found applications to a 
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FIG. 2: (color online) Rescaled propagators of Levy walk 
process with exponent 7 — 5/3 for different times, t — 
100, 300, 600. The dashed line depicts the propagator for the 
standard Levy walk with the constant flight velocity vo — 1. 
For the simulations of the random walk in the active medium 
we employed the map © with D„ = 0.03. The inset shows 
the ballistic front regions after the scaling transformation ((5]), 
where x — x — vot. 



wide range of phe nomena, ranging from financial mar- 
kets dynamics jl3j to single molecule spectroscopy 
Here, we use one specific CTRW model 0,111, where a 
walker moves ballistically in between successive "turning 
points" . During a single flight event, the walker trav- 
els at constant speed wq, and at the turning points it 
randomly changes the direction of motion (see the par- 
ticle flying above the slab on Fig. 1). Hence, the ve- 
locity probability density function (PDF) reads h{v) = 
[6{v — vq) + 6{v + vq)] /2. The flight times, t„, are in- 
dependent and identically distributed random variables 
drawn from a PDF ^{t) that is described by a power law 



tP{t) (X (t/to) 



(1) 



where tq is a characteristic time scale, and 1 < 7 < 2. 
This choice guarantees a flnite average flight time, (r) = 
Jp Tip(T)dT, and provides access to different diffusion 
regimes with the scaling of the mean squared displace- 
ment (J^{t)(x t^~^ . The corresponding Levy walk (LW) 
approach Q has been successfully used for the descrip- 
tion of diffusion of particles in chaotic systems 0, [13, [111 ' 
tracers in turbulent flows [l^l, or ultracold atoms in op- 
tical potentials 2l|. For the standard LW process, the 



PDF of finding a particle in x at the time t, P{x, t), pro- 
vided it was initially localized at x = 0, exhibits a sharp 
cutoff marked by the ballistic peaks at \x\ = v^t [l9| . 

A first step towards the implementation of the LW for- 
malism for the description of the perturbation dynamics 
in many-particle systems has been attempted in [Hj. Al- 
though this approach provided an adequate description 
for the perturbation spreading process, it failed to cap- 
ture the dynamics of the cone fronts. The observed ballis- 
tic humps displayed smooth, Gaussian-like profiles, with 
a scaling that was incompatible with the scaling behavior 
of deltalike peaks of the LW propagator 0. 



In order to resolve this issue, we here extend the con- 
ventional CTRW setup by assuming that the walker per- 
forms a random walk through an active medium. Con- 
ceptually, it means that while moving, the walker inter- 
acts with the surrounding medium (the slab in Fig. 1). 
This interaction causes fluctuations of the walker's ve- 
locity (note the noisy trajectory of the particle moving 
inside the slab in Fig. 1). The term "active" means 
that the medium is not solely dissipative and the parti- 
cle not only continuously loses but also gains energy from 
its environment. Both processes are in balance, yielding 
unbiased fluctuations of the walker's velocity around vq. 

Model dynamics. We start out by considering a single 
flight event. The position of the walker is deflned by a 
Langevin equation x = vq -\- where ^(i) is a delta- 
correlated Gaussian process of vanishing mean and finite 
intensity _D„, i.e., {£,{t)£,{s)) = Dy6{t — s). This consti- 
tutes a well-known biased Wiener process with drift vq 
[2^ . After an integration over a time interval r we obtain 



!(t + r) = x{t) + vqt + w{t), 



(2) 



where the new stochastic variable w{t) ~ J^'^'^ £,{s)ds 
is characterized by the Gaussian PDF p{w,t) with the 
dispersion cr^ = {[x{t) — wqt]^) = D„t. The propagator 
can be calculated numerically by sampling long enough 
a single-particle trajectory [23, (see Fig. 2). 

To gain analytical insight into the generalized LW dy- 
namics, we follow a standard reasoning [l|, Q and derive 
the transport equations governing the evolution of the 
particle density, P{x,t) (see supplemental material). Be- 
low we present the major results for the scaling properties 
of the central part of the density profile and provide the 
explicit expression describing the ballistic humps. The 
asymptotic analysis of the central part of the density pro- 
file reveals the scaling of the standard LW propagator [H , 
namely. 



P{x,t') 



1 



-P 



\x\<_vt (3) 



where K <x Tq '^^^va and u ~ t' /t, see Fig. 2. 

The salient difference between the dynamics of our 
model and the standard Levy walk becomes apparent 
in the regions of cone fronts. The ballistic humps are 
formed by the particles which were fiying from < = to 
the observation time t (see supplemental material): 

Phump(a;,t) = Ht)[p{x + VQt,t) + p{x - vot,t)] /2 (4) 

Here, $(t) denotes the probability of not changing the 
direction of flight during the time t and has a power-law 
behavior $(i) cx (t/ro)^"'*' 01 ■ Consequently, the area 
under the ballistic humps (jl]) also scales as t^~'^ . During 
flights, the particles undergo random fluctuations caused 
by the flight's velocity variations. The flight length is 
proportional to t; thus the dispersion of the Gaussian- 
like humps grows as y/t^ and we arrive at the following 
scaling for the particles' density in the hump regions: 

Phump(5, t') ~ U-^/^P^^^^ix/u-^-^/^ t), (5) 
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and periodic boundary conditions. Without loss of gen- 
erality we set the mass ratio M/m = 2. The energy per 



FIG. 3: (color online) (a) Rescaled perturbation profiles at 
times t = 1000, 2000, 4000, and 6000 (the width increases with 
time), for the hard-point gas dynamics with the energy per 
particle e = 1. The scaling exponent is 7 = 5/3. The inset, 
similar to the inset in Fig. 2, depicts the ballistic humps after 
the scaling transformation ([Sjl. (b) Profiles after the scaling 
transformation (jSj for different energy per particle, e, at time 
t = 1500/ i/e/eo, with eo = 1- Each profile is obtained by 
averaging over 10^ realizations. 



where u = t' /t and x = x — vot (see inset in Fig. 2). Note 
that this scaling distinctly differs from the scaling in Eq. 

m- 

Ergodic many-body systems: validation of the ap- 
proach. Consider a many-particle system, with a Hamil- 
tonian 



N 



Htota.l{,{Xi,Pi}) = y^iJi 



(6) 



where Hi = H{xi,Xi-i,Xi+i,pi) is the energy attributed 
to the ith particle. At the time t = the system is lo- 
cally affected by the perturbation. The initially localized 
perturbation energy, Ep, starts to spread, such that the 
distribution of the local excess energy AE{i,t) of the i 
particle evolves in time [25l |. while keeping the pertur- 
bation energy constant; i.e., E^j^Ai?(i,t) = Ep. The 
spreading can be quantified with a normalized probabil- 
ity distribution function g{i,t) = AE{i,t)/Ep, where 
denotes a microcanonical average. 

The main finding of this study is that the profiles of 
the spreading perturbation in many-particle Hamiltonian 
systems, g{i,t), in the corresponding asymptotic regime 
|26| . are determined by the propagator P{x, t) of the gen- 
eralized Levy walk model, Eqs. (fT2l [5|). In order to val- 
idate this claim we use two archetype systems, namely 
(i) a one-dimensional scattering dynamics of a hard-point 
gas composed of two unequal masses [27^ an d (ii) a Fermi- 
Pasta- Ulam (FPU) /3-lattice dynamics [28 1. 

We start with a hard-point gas, a many-body Hamil- 
tonian system with an ergodic dynamics governed by the 
conservation of kinetic energy and momentum (27| . We 
use a chain of = 1.6 • 10^ pointlikc particles with alter- 
nating masses, ...mMm,M..., of the length L = 1.6 • 10"*, 



particle, e 



m or M , serves as a tun- 



able parameter. Figure 3 depicts the evolution of the 
infinitesimal perturbation Q{i,t) Q. The scaling ansatz 
([5]) with the exponent 7 = 5/3 is beautifully validated, 
see the inset in Fig. 3a. We also found that the pertur- 
bation profiles for different values of energy per particle 
parameter e are matched by assuming that the pertur- 
bation velocity and the fluctuation variance both scale 
as 



Consequently, the profile scales as 

Qe{x,t) = ge'{x,t/s'), 



(7) 
(8) 



where s' = \/e'/e, see Fig. 3b. 

As our second test bed we use a FPU (3 chain dy- 
namics [2^, defined by the Hamiltonian ^ with Hi = 
iPi + \{xi+i - Xi)"^ -l-f (xi+i - XiY, with N particles 
of unit mass and periodic boundary conditions. The en- 
ergy per particle is e = iJtotai/A^- It is not feasible to 
explore the evolution of finite perturbations of the FPU 
system at microcanonical equilibrium, due to emerging 
huge statistical fluctuations. Instead we employed the 
energy correlation function, e{i,t) @, 0], which bears the 
same information as the infinitesimal perturbation in the 
case of hard-point gas 29| . We performed a massive nu- 
merical experiment [30 1, yet even these efforts were not 
sufficient to cope with the statistical fluctuations [note 
the thin green lines in Figs. 4a, b]. Nevertheless, a rel- 
atively smooth shape of the ballistic hump allows for a 
convincing validation of the scaling . It is interesting 
to note that the velocity of the ballistic peaks is deter- 
mined by the group velocity of effective thermal phonons 
0, which, therefore can be associated with "walkers" of 
the CTRW approach. 

The perturbation profiles for different values of the en- 
ergy per particle e reveal another remarkable feature: the 
central part of profiles is independent of e, while the bal- 
listic humps move as e increases, see Fig. 4c. Such be- 
havior can be derived from the scaling invariance (|12l) . 

Taking into account that if oc Tq ^^^uq, one can infer 
that the central part of the LW propagator is invariant 
under variation of wq when tq oc tij 

In conclusion, we demonstrated that the collective 
process of perturbation spreading across two celebrated 
many-particle systems, a hard-point gas and a Fermi- 
Pasta-Ulam lattice, is reproduced by a single-particle 
stochastic process. It is intuitive that an ergodic dy- 
namics is a prerequisite for the diffusionlike perturba- 
tion evolution. In integrable many-bod y s ystems, such 
as a harmonic chain or the Toda lattice |3l| . the pertur- 
bation spreading is a deterministic process, evolving in 
terms of nonintcracting phonons. The task of exploring 
the sufficient conditions for the CTRW kinetics to occur 
presents a promising challenge, thereby underpinning the 
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FIG. 4: (color online) (a, b) Energy correlation functions 
e{i,t) for the FPU dynamics with iV = 1.6 • 10*, ^ = 1 and 
e = 1 (thin solid lines), compared to the propagators of the 
generalized LW model with vo = 1.384 and Dv = 0.49 (thick 
dashed lines). The scaling exponent is 7 = 5/3. The FPU 
profiles have been obtained by averaging over 5 ■ 10"" realiza- 
tions, (c) Energy correlation functions for the FPU system 
with N = 2 ■ 10"^ for the time t = 500 and different values of 
the energy per particle e. 



universality of our findings. Thus our results disclose a 
pathway to explore propagation of information in real- 
istic dissipation-free systems: it allows to calibrate the 
transport characteristics of many-body systems (which 
are beyond the region of validity of the classical Lieb- 
Robinson theory [ll|) by using the parameters of the 
corresponding random walk process. 



Apart from those theoretical challenges, there is room 
for possible applications. A realization that comes to 
mind is an array of coupled nanoresonators js^ . where 
a single unit acts as both, the receiver and transducer of 
excitations, which transforms this array into an extended 
sensor via utilizing the principle of time-of-arrival local- 
ization [33j . 
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I. TRANSPORT EQUATIONS FOR THE 
MODEL PROPAGATOR 

In order to gain analytical insight into the general- 
ized LW dynamics, we proceed along known reasoning, 
see in [l|, and derive the transport equation for the 



probability distribution function (PDF) P{x,t). We first 
introduce the probability distribution of the end points 
of flights, or the turning points, i>{x,t), where a particle 
chooses its new velocity. Its evolution is governed by the 
following balance equation: 



i/(a;, t) 



dvdw / v{x ^ VT — w^t — T)il){T)h{v)p{w,T)dT -\- tp(t) I h{v)p{x — vt,t)dz 



(9) 



r 



Equation ^ shows that a particle changes its velocity 
at the point {x,t) if it was the end point of the pre- 
ceding step. That previous step had some flight time r 
and occurred with some velocity v. Therefore the step 
originated in the point x — vt — w, where w takes into 
account the accumulation of velocity fluctuations dur- 
ing a single flight. Next we have to integrate over all 
possible flight times, velocities and fluctuations with the 
corresponding probability densities. This is how the first 
term on the right side of Eq. ^ is obtained. The second 
term reflects the influence of initial conditions. In this 
work we use so-called equilibrated initial conditions 



which assumes that walkers evolved for an infinitely long 
time when the observation started. Different initial se- 
tups affect only the probability of when a flying particle 
experiences the first turn after the start of observation. 
For a system evolving for an infinite time the PDF of the 
first turn after the observation has been started is given 
by Lp{t) = (t)~^ Jq°° ilj{t + T)dT 0]. If a particle starts at 
X = 0, the spatial position of where the first turn occurs 
is infiuenced by the velocity fluctuations and is given by: 
loo — vt — w)p{w, t)dw = p{x — vt, t). 

We next evaluate the actual density of particles P(a;, t), 
to obtain: 



P{x,t) = J J dvdwjv{x-VT-w,t-T)'^[T)h{v)p{w,T)dT + '^{t) / h{v)p{x - vt,t)dv . (10) 

— OO — OO 

I 



The particle currently flnds itself in the point {x,t), if 
it previously changed the direction of flight at the point 
X — VT — w, and keeps flying afterwards for the time r 
with the probability ^'(t) = 1 — tjj{t)dt. Similarly, 
the probability to continue the flrst flight reads $(i) = 
(T)-Vo°°*(t + r)dr i. 

Formally, the above two equations can be solved 
with a help of combined Fourier- and Laplace-transform, 
thereby turning all convolution integrals into products 
and thus rendering the integral equations algebraic. The 
answer for the density of particles in Fourier/Laplace 
space, Pk^s is given by (we use tilde-notation to denote 
the Fourier/Laplace transform): 



*(T)ft.ferPfc(T") ViT)hkrPk{T) 



1 - [/lfcrPfc(T)-0(T)] 



Pk,s = : -+ hktPkirMt) 

(11) 

This exact analytical expression serves as the starting 
point for the asymptotic analysis for large spatial and 



temporal scales, that corresponds to small fc, s coordi- 
nates in Fouricr/Laplacc-spacc. It is possible to show 
that for the case of Gaussian fluctuations the central part 
of the proflle can be described by the Levy distribution, 
which is the solution of the fractional diffusion equation 
Hi: 



Pk 



S + T-n' 



v^{j - l)r[l - j]k'i sin(7r7/2) 



In original coordinate space and original time this expres- 
sion delivers the following scaling relation for the central 
part of the density proflle: 



P(x,i') 



1 



'(^,^), \^\«vt, (12) 



where K <x ^^''^wq and u ~ t' /t, see Fig. 2. 

To describe the density of particles in the ballistic 
humps of the profile we use the second term on the right 
hand side of Eq. pi!]) : 
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II. NUMERICS 



Phump(a;, t) - $(0 J h{v)p{x - vt, t)dv , (13) p^^. ^j^g integration of the FPU /3 chain's equations of 

-oo motion we used the symplectic SABA2C scheme [H, 

with the integration time step dt = 0.01 ~ 0.02. Calcu- 
where $(<) = (t)~ + r)dr. For the chosen flight lations have been performed on a Tcsla S1070 supercom- 

time distribution it scales as <^{t) oc (i/ro)^"'''. puter, with 960 CPU's on board. 
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